Stochastic differential equations for non-linear hydrodynamics 
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We formulate the stochastic differential equations for non- 
linear hydro dynamic fluctuations. The equations incorporate 
the random forces through a random stress tensor and ran- 
dom heat flux as in the Landau and Lifshitz theory. However, 
the equations are non-linear and the random forces are non- 
Gaussian. We provide explicit expressions for these random 
quantities in terms of the well-deflned increments of the Wi- 
enner process. 



I. INTRODUCTION 

Two physical phenomena affect in a crucial manner 
the dynamics of concentrated submicron colloidal sus- 
pensions. On one hand, the hydrodynamic interaction 
between the colloidal particles is responsible for the com- 
plex behavior of these systems. The motion of a colloidal 
particle affects the state of motion of the other particles 
through the perturbation of the mediating solvent. Be- 
cause of the long range nature of these interactions, this 
is essentially a many-body problem for which analytical 
results are very difficult to obtain. The second physical 
phenomena that intervenes crucially in the dynamics of 
these systems is Brownian diffusion that keeps the parti- 
cles in constant thermal agitation. 

Hydrodynamic interactions and Brownian diffusion are 
the macroscopic and mesoscopic manifestations of the hy- 
drodynamic behavior of the solvent in which the colloidal 
particles are immersed. Actually, the origin of the Brow- 
nian diffusion can be traced back to the fluctuations of 
the hydrodynamic fields surrounding the particles. This 
explanation of Brownian diffusion constitutes the hydro- 
dynamic theory of Brownian motion which has 
been successfully used to predict the celebrated long-time 
tails in the velocity autocorrelation of a single Brownian 
particle. This prediction has been confirmed by experi- 
ments @- @. 

For the case of undiluted suspensions there is no ana- 
lytic theory that describes simultaneously the rheological 
and diffusive aspects of the system. It seems that the 
only viable way is to recourse to numerical simulations. 
For this reason, there is at present great interest in meso- 
scopic levels of description for complex fluids that allow 
to address such problems with affordable computation 
work. In this respect, lattice Boltzmann automata 
I p^ and dissipative particle dynamics |18| have been 
useful tools in studying hydrodynamic problems as those 



occuring in immiscible fluids and colloidal and polymeric 
suspensions. 

In the lattice Boltzmann technique for colloidal sus- 
pensions, one solves a simplified version of the fiuctuating 
Boltzmann equation on a lattice The advantage of 
the technique is that, due to the local description of the 
hydrodynamic behavior, the computational time grows 
linearly with the number of colloidal particles instead of 
cubically as in Stokesian dynamics |l9| . The drawback 
lies in the difficulty of handling irregular moving bound- 
ary conditions on a lattice. On the other hand, the spirit 
of dissipative particle dynamics is to formulate point par- 
ticles that move similarly as the "droplets" of the fluid 
would move, thus capturing the hydrodynamic behavior 
in an off-lattice technique. This approach is similar to 
the smoothed particle hydrodynamics technique [pO 21 
in which the macroscopic hydrodynamic equations are 
discretized on moving points. Smoothed particle dynam- 
ics, as opposed to dissipative particle dynamics, does not 
include any thermal randomness into account. There- 
fore, it is not clear how the fluctuations induced by the 
discreteness of the system are consistent with thermal 
fluctuations that have a microscopic origin. In order to 
clarify this point it is desirable to have a well-defined for- 
mulation for the Langevin equations for hydrodynamic 
fields. 

There exists a well-known theory for hydrodynamic 
fluctuations which was formulated by Landau and Lif- 
shitz in 1957 ||2^]. It consists basically on adding stochas- 
tic fluxes to the stress tensor and heat flux. The noise am- 
plitude is determined by the temperature and the trans- 
port coefficients of the fluid. The main problem in try- 
ing to apply the Landau-Lifshitz theory for numerical 
simulations of time-dependent flows as those arising in 
fiuctuating fluid solvents in colloidal suspensions can be 
summarized in the question, what value of the tempera- 
ture (fluctuating and in general non-uniform) should one 
use in the noise amplitude? And, if the transport coef- 
flcients depend on the state of the fluid, what values of 
the hydrodynamic variables should one use, fluctuating 
or average values? These questions should be answered 
satisfactorily in any attempt to conduct numerical sim- 
ulations of non-equilibrium time-dependent fluctuating 
fluids. 

The Landau and Lifshitz theory deals with linear fluc- 
tuations around equilibrium states but the theory has 
been further generalized to linear fluctuations around 
non-equilibrium steady states However, hydrody- 

namics contains non-linearities coming from the depen- 
dence of the transport coefficients on the state variables 
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and to the presence of quadratic terms in the Raighley 
dissipation function. There have been several attempts 
to study non-linear hydrodynamic fluctuations within 
the framework of Fokker-Planck descriptions j^J]- [^ . 
Fokker-Planck equations (FPE's) for hydrodynamic fluc- 
tuations are well-known since the early 80's. What does 
not seem to be known are the stochastic differential equa- 
tions (SDE) associated to the FPE. In the most general 
and complete work dealing with FPE for hydrodynamic 
variables Zubarev and Morozov conclude that "the prob- 
lem of the adequate Langevin equations describing non- 
linear hydrodynamic fluctuations seems to be still open" 
p8| . In another important paper, van Saarloos, Bedeaux 
and Mazur formulate the Langevin equations for non- 
linear fluctuations by assuming that the random fluxes 
are Gaussian processes with white noise |27 . This leads 



to an inconsistency in the form of particular unrealistic 
functional dependences of the transport coefficients on 
the temperature. In addition, the Einstein distribution 
function is not the equilibrium solution of the correspond- 
ing Fokker-Planck equation . 

In this paper we formulate the appropriate SDE for 
nonlinear hydrodynamic interactions starting from first 
principles, i.e. from the FPE for hydrodynamic fields 
derived in Ref. The derivation of these Langevin 

equations is a necessary initial step in order to have well- 
defined starting points for numerical algorithms in which 
both the macroscopic and mesoscopic aspects of hydro- 
dynamics are captured consistently. We provide explicit 
expressions for the random stress tensor and the ran- 
dom heat flux in terms of the increments of the Wienner 
process. The fluxes are shown to be non-Gaussian pro- 
cesses. These explicit forms for the random fluxes allow 
to implement a numerical simulation of the fluctuating 
hydrodynamic equations straightforwardly. 

The formulation of the stochastic differential equations 
as appear in this paper has an additional advantage when 
performing numerical simulations over the usual linear 
hydrodynamic fluctuations approach |Q. For the 
case of linear hydrodynamic fluctuations around non- 
equilibrium states one has to first compute the non- 
equilibrium state of the system (usually assumed to be 
stationary) and then linearize the equations around this 
steady state. This involves a two-step procedure whereas 
in the formulation of the present paper a single simula- 
tion produces simultaneously the mean motion and the 
fluctuations in one go. This is similar in spirit to the 
hydrostochastic formulation of Breuer and Petruccione, 
where a master equation is formulated for the hydrody- 
namic fields [^ . It is apparent that for the simulation 
of the hydrodynamic solvent between colloidal particles 
where the non-equilibrium (non-stationary) state is not 
even known, the approach of linear hydrodynamics is 
completely useless. 

A word is in order about the procedure we follow for 
deriving the stochastic differential equations for the hy- 
drodynamic equations. It is possible to obtain from first 
principles (for example with projection operator tech- 



niques) the Langevin equations governing the hydrody- 
namic variables. However, we have discussed at length 
in Ref. that this procedure has the shortcoming that 
it does not provide information about moments of the 
random forces higher than the second one. Therefore, 
issues about the Gaussian nature of the noise cannot be 
resolved. For this reason, we follow in this paper a differ- 
ent route that is based on the mathematical equivalence 
between the Fokker-Planck equation and the correspond- 
ing stochastic differential equation . Given a Fokker- 
Planck equation one can read from the coefficients of the 
equation which is the corresponding Langevin equation 
for the stochastic process (with its appropriate stochas- 
tic interpretation, Ito or Stratonovich). There are no 
physical approximations involved in this procedure. A 
close analogy arises with the connection between the Li- 
ouville equation and the Hamilton equations in classical 
mechanics. They are mathematically equivalent repre- 
sentations of the same process (in this case a determinis- 
tic process). The actual interest of looking at Langevin 
equations instead of Fokker-Planck equations arises from 
a computational point of view. Langevin equations are 
much more easier to simulate than to solve a partial dif- 
ferential equation as the Fokker-Planck equation. 

This paper is notationally a bit cumbersome. This is 
quite an inescapable situation given that we are dealing 
with five continuous hydrodynamic fields. In order to 
alleviate a bit this we have decided to give the struc- 
ture of the Fokker-Planck equation and the correspond- 
ing stochastic differential equation for arbitrary fields be- 
fore particularizing to the case of hydrodynamic fields. In 
any case, the final results have a very compact structure 
closely reminiscent to the Landau and Lifshitz theory. 

The paper is distributed as follows. In section II we 
discuss the definition of the coarse-grained hydrodynamic 
variables in some detail. In section III we rewrite the 
standard FPE for the case that the set of relevant vari- 
ables are extensive variables characterized by a discrete 
index associated to a lattice (i.e. a field defined on a 
lattice). We also write the FPE equation corresponding 
to the densities of the extensive variables. This FPE is 
then expressed in terms of functional derivatives. Sec- 
tion IV is devoted to the equivalent SDE for fields on a 
lattice and its corresponding functional Langevin equa- 
tion where precise meaning to the white noise term is 
given. In section V, the general FPE for fields on a lat- 
tice obtained in section III is particularized to the case 
that the fields are the hydrodynamic fields. Section VI is 
devoted to the derivation of the SDE for hydrodynamic 
fields. Finally, we discuss the results in the last section. 



II. HYDRODYNAMIC VARIABLES 

Our aim is to obtain the FPE for the hydrodynamic 
continuous fields. From a mathematical point of view 
there are subtleties regarding the existence of the con- 
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tinuum limit for stochastic variables and for this reason 
we will restrict ourselves to the case that the fields are 
defined on a lattice. In this section we discuss the dis- 
crete variables that are representative of the hydrody- 
namic fields. A word is in order about the definition of 
the lattice of points r. This lattice is usually taken as 
a regular cubic lattice ||2^. In the formulation that fol- 
lows there is no restriction, however, to the kind of lattice 
that can be used and irregular lattices are also possible. 
To each lattice point we associate a volume element of 
volume Wr that defines the region of space associated to 
each lattice point. In the case of the cubic lattice it is 
natural to associate a cubic cell of volume e'^. For ir- 
regular lattices, the Voronoi tessellation allows to define 
non-overlapping volume elements that cover all the space 
. Note that the cubic cell is a special case of Voronoi 
tessellation for the case of a regular cubic lattice. 

We take as relevant variables the total mass Mr, the 
total momentum Pr and total energy of that set of 
particles that are within the volume element associated to 
the lattice point r. In compact form, Mr,Pr,£^r — ^ ^™ 
with m — 0, ... ,4. Because these variables are conserved 
variables, one expects that their dynamics is much slower 
than any other set of dynamic variables in such a way 
that they form a closed description for the system. The 
size of the volume elements should be so large that the 
relevant variables can be considered as continuous. For 
example, if a single particle enters the cell, the variation 
of the mass of the cell can be considered as infinitesimal. 

In mathematical terms, the relevant variables are given 
by 



Mr ^^mXr{(l^) 

i 

i 



(1) 



where we denote by qi,Pi the position and momentum 
of particle i and its energy is 



Pi 

2m 
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(2) 



The function Xr(q) is the characteristic function of the 
cell centered at r, which takes values 1 if q belongs to 
the cell and otherwise. The volume of the cell at r is 
given by 



dr'xr(r') 



(3) 



If we displace the position q^ of all the particles a vector 
r' and simultaneously displace the location of the cell at 
r the same vector r', the numerical value of Xr(qi) will 
not change. This translational invariance implies that 



Xr(qi) = Xr(qi — r). It is apparent that if all cells are 
identical in shape as occurs in regular lattices we can 
suppress the label r of the characteristic function, i.e. 

Xr(qj - r) = x(qj - 1"). 

The density variables associated to our selected exten- 
sive variables are 



AL 



= ^mSr{r - qi) 



gr = — = Vpj5r(r- qi) 

-J 7.- ^ 



Cr = — = BiSrir - qO 



(4) 



where we have introduced 

Mr') = Xr{r')/vr (5) 
that satisfies by virtue of the normalization condition 

dr%{r') = 1 (6) 



(7) 



The velocity field is defined by 



Note that X^r ^rPr ~ Sr ^'^r, stc. Wc want to have 
J2r -^^r squal to the total mass of the systems and this 
is only possible if the cells are non-overlapping and cover 
all the space. If this condition is not met, then Mr 
is not a conserved variable and one does not expect it to 
be a slow variable. 

For the case of a cubic lattice with cubic volume ele- 
ments, the characteristic function can be represented as a 
product of square step functions in the x,y, z directions. 
One can argue that if the particles have some finite size, 
then it is meaningful to smooth out the step functions at 
the borders of the cell and to talk about fractions of a 
particle being in different cells. Actually, one could use 
coarsening functions other than i5r(r') or these smoothed 
square step functions defining the "shape" of the cell. 
Zubarev and Morozov make use of a particular coarse- 
graining procedure in which the hydrodynamic fields are 
represented in a Fourier basis and then only the small 
k < ko components are retained . In real space this is 
equivalent to use a coarse graining function of the form 



X{r) 



l3 

27r2 



sin fegr cos fcnr 



■0' 



(fcor)3 (fcor)2 



(8) 



This particular coarse-graining produces a coarse-grained 
density field that is not everywhere positive. This can be 
cured, for example, by using a Gaussian function 



A(r) 



1 



7r3/2 



exp — (r/e) 



(9) 
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where e is a coarse-graining length. Nevertheless, in this 
paper we will restrict ourselves to the use of the coarse- 
graining function (5r(r') derived from the characteristic 
function of each cell, Eqn. (^). 

We will need the time derivatives of the hydrodynamic 
variables. They have the form 



(10) 



where the coarse-grained stress tensor and energy flux 
are given by the usual expressions (see, for example, Ref. 

0)- 

A remark on the meaning of Vr is in order. In obtain- 
ing (ll^) we have used 



4^r(r-qO = -|:^r(r-q,) 



(11) 



where we assume that the index r is continuous. In fact, 
Eqn. ( [TTI ) is a mathematical identity. However, we will 
evaluate the right hand side of (|l^) only at the discrete 
values of the vector r. There is no need to use discretized 
forms for the nabla operator with this understanding P7| . 



III. FPE FOR FIELDS 

Zwanzig derived from first principles the Fokker- 
Planck equation for the distribution function of a discrete 
set of variables with the aid of a projection operator tech- 
nique |]35|| . The derivation provided explicit expressions 
for the drift and diffusion terms of the equation in terms 
of microscopic quantities. In this section, we will present 
the FPE for the case that the set of discrete variables are 
extensive variables A™ (or simply A) defined in each cell. 
We will also give the FPE for the distribution function 
of the densities = /vr- 

Let us consider an isolated system governed by Hamil- 
ton's equation of motion {H{z) is the Haniiltonian and 
iL is Liouville's operator for the dynamics). The micro- 
scopic degrees of freedom are denoted by z. The FPE is 
a differential equation for the probability density p{a^ t) 
that the system takes values a at time t for the phase 
hmctions j4(z), i.e. 



p{a,t) = / pt{z)S{A{z) — a)dz 



(12) 



where 6{A{z) — a) is a product of Dirac's delta functions, 
one for each relevant variable, and ptiz) is the solution 
of Liouville's equation. 
The FPE is [ESlJslaslJSTtl 



dtpia, t) 



E 



d 



2 E E 

mr m'r' 



d d 



with the following definitions 

K^\a) ^ <(a) + ^ C"'(«)^"'(") 



E 



dal 



/•mm / \ 



(14) 



with 



C"'(«) 



iLA';Y 

du{{SiLA'^5iLA^'{u))°' 



d 



da" 



lnr2(a) 



5{A{z) — a)dz 



/ ^ \ /-mm \ I /-m m / „ \ 

(a) = Crr' (a) + Cr'r (") 



(15) 



where we have defined the projected part SiLA^ = 
iLA^ - (iLA^)". 

We have introduced the average of an arbitrary phase 
function G over the submanifold A{z) — a by 



1 



S{A{z) - a)G{z)dz 



(16) 



The derivation of the FPE from the underlying micro- 
scopic dynamics of the system relies on the fact that 
the relevant variables A are much slower than the rest 
of variables of the system. This allows for a Markovian 
memoryless description of the dynamics of these relevant 
variables. A systematic expansion of a microscopically 
exact master equation in terms of a slowness parameter 
(which is, essentially, the ratio between time scales be- 
tween relevant and irrelevant variables) has been carried 
out by Mori et al. To second order in the slowness 
parameter the non-linear FPE ( |l3| ) is obtained. 

In the microscopic formulation of the FPE the drift 
and the diffusion tensor are defined in terms of 

microscopic quantities. The problem is that these quan- 
tities are difficult to compute. Basically, three difficul- 
ties have to be resolved: the calculation of the averages 
(• • •)", the calculation of the thermodynamic forces F^ 
(or equivalently of the phase space volume ^{a)) and fi- 
nally the calculation of the correlation functions in the 
kinetic tensor C™,™ . The first two are closely related and 
can be resolved for the case of hydrodynamic variables by 
resorting to the local equilibrium approximation. The ex- 
plicit calculation of the kinetic coefficients is much more 
difficult and they are simply left as phenomenological co- 
efficients. 

The rest of the section is devoted to the FPE equa- 
tion for the distribution function of the densities A™ = 
^rV^r- The corresponding FPE is obtained by a simple 
change of variables 



(«M«,i) (13) dtp{Vr&,t) 



mr 



d 1 
da™ Vr 



K^{vra)p{vra, t) 



4 



1 



d d 1 



mr m'r' r i 



(wra)p(wrQ;,t) 
(17) 



where a™ = /vr- We therefore introduce in a natural 
way the following quantities 



M 



p[a,t 



W_Vrp(va,t) 



—K^{va) 

Vr 

n{va) 



(18) 



We use square brackets to denote that the left hand sides 
have different functional forms from the right hand side 
(for example, is a different function of its argument 

from F^{va) even though they are named with the same 
symbol F™). M is the total number of variables. In this 
way the probability density P[a,t] is normalized to unity, 
i.e., J VaPla, t] = 1, where Va is a short hand for da^^ . 

Note that these definitions (|l^) allow to express the 
reversible part of the drift and the kinetic tensor in terms 
of averages over density fields 



du {{6iLA'^){6iLA"J {u))Y^ (19) 



The FPE for the densities ( poj ) admits a suggestive con- 
tinuum notation. By making the substitutions of partial 
derivatives by functional derivatives according to 



1 d 

Vr 9(5™ 

and sums by integrals 



5cC 



dr 



we can write equation ( ^ ) in functional form 
dtP[&,t] - - Jdr^K^[&]P[a,t] 



with 



Sa") 



T^mm r~i /•mm r~i , /-mmr~i 
Drr' m = Crr' ["J + Cr'r W 



(22) 



(23) 



(25) 



The functional FPE ( P4[) is understood in this work sim- 
ply as a notational device and the mathematically well- 
defined FPE is the discrete FPE (|^. 



IV. STOCHASTIC DIFFERENTIAL EQUATION 
FOR FIELDS 



where we have defined the projected part diLA™ = 
iLA^' - {iLApY^. With the definitions (|l|) the FPE 
for the density fields is 



dtP[a,t] 



with 



E 



d 



■x;"[(5]p[(5,t] 



\llY.iLii^D-r\^]n^A (20) 



mr m'r' 



7^, 9(5™ 9(5™ 



r' 

9 



as™' 

1 9 



■lnl7f(5l 



Wr 9(5™ 

Tiram r-l /-rnm r~i i 
-Drr' ["J = Crr' ["J + C 



(21) 



As it is well known associated to the FPE there 
is a stochastic differential equation (SDE) describing the 
stochastic process a{t). The SDE corresponding to the 
FPE is given by 



d(5;"(t) = if;"[(5(t)]di 



mm [~r 



a(t)]dVKr7 (i) (26) 



where we have introduced the independent increments of 
the Wienner process associated to each cell r. These in- 
crements satisfy the mnemotechnical rule for Ito calculus 



dW^\t)dW^ {t') = 6^,n'Srr'dt 



if t = i' (27) 



and zero ii t ^ t' . The Kronecker delta 5rr' is one if 
r = r' and zero otherwise. The matrix _B™™ [5] is the 
square root in matrix sense of 13™™ [5] , this is 



(28) 
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The SDE (^6|) is interpreted in Ito sense. This inter- 
pretation is most convenient for numerical simulations 
because the increment of the Wienner process is uncor- 
related with the variables at the same time. For a discus- 
sion of the advantages of deriving stochastic differential 
equations from first principles through the FPE route in- 
stead of directly from the microscopic dynamics, see Ref. 

A physically suggestive way of writing the SDE (g6|) is 
by using white noise in space and time defined as 



77'"(r,i) 



(29) 



We thus write Eqn. (|2^) as an ordinary differential equa- 
tion, i.e. in the form of a Langevin equation 



dta'^it) ^ K^[a{t)] + ydr'B™r [5(t)]?7™ (r',t) (30) 

where we have introduced the "continuous" matrix 
B^V'[ait)] = B™,™'[a(t)]/7i^ that satisfies 



(31) 



The white noise term has the following correlation 
dW^{t) dW^'it') 



(rr(r,t)77'"'(r',0> = ( 



= Sn 



Vj-dt ,Jv^dt' 
Vj. dt 



iit = t' (32) 



and zero if t =/= t'. We can rewrite this correlation as 
follows 

(77™(r,i)77'"'(r',t')) ='5,w5(r-r')5(t-i') (33) 

The reason for such a procedure is that expression (^) 
will give the same results as (^2|) when computing aver- 
ages and correlations. Indeed, let us consider the noise 
term in (Eq) which is given by 



/™(t)dt eeJ^B^v^' [a{t)]dW"''ir,t) 



(34) 



The notation suggest that /"(t) is related to the deriva- 
tive of the Wienner process that does not exist in a strict 
sense. Nevertheless, with the understanding that f^{t)dt 
will appear always inside integral expressions (which will 
be interpreted as stochastic integrals) such a procedure 
leads to correct results. For example, the correlation of 
the noise terms (Bif) will be 



dt dt"r:\t)ff{t") 



E 



i(i)]dM^"'(r',i) 



X B'p,"X' [a{t')]dW"'"' (r'", t')j 
= [ dt [ VaP[a,t]D^;?"[a] = / dt/u^'JV") (35) 



where Ito rule (|3) has been used. Note that one can 
write 



/r(t)/;s it") ) = 6{t- 1") (i^™" )^ (36) 



because if one integrates over t, t" one recovers ( p5| ) and 
the integration limits are arbitrary. 

On the other hand, in the continuous version ( ^ ) the 
noise term is given by 



/r(t)^ / dv'B'^^r [<5(t)]77'"(r',t) 



(37) 



Now, if we use the continuum version ( |37| ) and proceed 
with the ordinary calculus rules (and the fact that the 
noise is uncorrelated with the variables at the same time) 
we will have 

ut Jdt'f^{t)f^:'it')^ 

(dt Jdt' J dr'B™r'[a(t)]r/™'(r',i) 
X / dr'''B™r"''''['5(i')]?7"''' (!■''', i') 

Idt Jdt' J dr' J dr"'(B™"'[5(t)]B™;™"'[5(i')] 
X (?7"'(r',t)r/™"'(r"',t')^ 

fdt I dr'(B™"'[<5(0]SX"'[fi(i)]; 

^ Jdt(^D^;/^'[a{t)]) (38) 

which provides the same result as (p5|). 



V. FOKKER-PLANCK EQUATION FOR A 
SIMPLE FLUID 

Now we wish to specify the FPE (|2^) for the case that 
the continuum variables are the hydrodynamic vari- 
ables. This has been done explicitly in Ref. ||2^ for the 
case that the coarse-grained hydrodynamic fields are de- 
fined in terms of the small k < ko Fourier components. 
In order to derive the SDE ( |30| ) we found more con- 
venient to work with the coarse-grained hydrodynamic 
fields defined in terms of the densities (i.e. extensive 
variables divided by the volume of the cell). For this 
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case the derivation of the exphcit forms for the drift and 
diffusion terms foUows very similar Unes as those of Ref. 
psf and therefore, only the final results are quoted here. 
The physical assumptions underlying the derivation of 
the Fokker-Planck equation for hydrodynamic variables 
are essentially that the fields are slowly varying in space 
and time. Therefore, a second order expansion in gradi- 
ents is allowed along with a Markovian assumption. Con- 
sistently, the concept of local equilibrium is advocated. 
These assumptions are thoroughly discussed by Zubarev 
and Morozov 

The thermodynamic forces i^™[a] are given by 



F;:[a] 



1 

T'r 
r 



(39) 



where is the temperature and /ir is the chemical poten- 
tial which are given in terms of the equilibrium equations 
of state, this is 



(40) 



Note that the 



where the internal energy is e = e 
thermodynamic parameters i^™[Q;] which are, in princi- 
ple, functionals of the fields become local functions of the 
density variables, this is, ^'^'"[q^] = ^""(ar)- 

The reversible part wj? [a] of the drift term ( p^ is given 

by 



(iier)" = -V-[Vr(er+Pr)] 



(41) 



where is the equilibrium pressure evaluated at the 
instantaneous values of the hydrodynamic fields, this is 



The kinetic tensor is given by (|19| ). Note that the com- 
ponents corresponding to the density field vanish because 
this field has no projected part, i.e. {iLp^)" = iLp^, see 
( p] ) . This imply that the only components different from 
zero correspond to the momentum density and the energy 
density. The explicit forms are 



Cr1=9r,ar'/(r-r')L?^'3''v^ 



(42) 



Lf^'' = T(r)77(r)5„^^, + T(r)C(r),S„;3^, 

...... 



(43) 



where ?7(r),C(r) are the shear and bulk viscosities and 
K(r) is the thermal conductivity. These transport coef- 
ficients are given in terms of the usual Green-Kubo for- 
mulae where the equilibrium averages are evaluated at 
the instantaneous values of the temperature and chemi- 
cal potential fields. 

The drift term if™ [a] in ( |25| ) is obtained by collecting 
the results for the reversible part w™[a] in (|l|), for the 
kinetic coefficient (^ ) and for the thermodynamic forces 
). We will have 



VrVrPr + Prl - 2777 ^ (C " o'?)tr[7]l 



Kr[a] = —dr- 



Vr(er 



Pr) - 277Vr-7 - (C - g?7)tr[7]Vr 



(44) 



where 7=5 [Vv + (Vv)"^] is the strain rate tensor. 

It is shown in Ref. Q that the term of K^[a] in (|2|) 
vanish for the case of hydrodynamic variables. This is 



_6_ 

56C' 



dr' 7 Cri-' 

c ~,m' ^rr 



mm r ~ 



[<5]=0 



(45) 



This is a consequence of the locality in space assumption 
and has the very important consequence that the Ito and 
Stratonovich forms of the FPE coincide . In this way, 
the resulting SDE can be interpreted in both senses with 
identical results. 

The diffusion matrix Z?™™ (p5|) is given by 



rr' 



r 

a^f5u 3 



where we have introduced 

+ T{v) Q,7(r)+C(r)) [S^.S^^ + S 



(46) 



(47) 



where we have introduced 



which is a tensor symmetric in a, /3 and pLjV. 

Now that we have the drift term and the diffu- 

sion matrix Z?™7' [a] we could write down the FPE (|20| ) 
explicitly. 
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VI. SDE FOR SIMPLE FLUIDS 

The main difficulty in writing down the SDE associated 
to the FPE for hydrodynamic fields is that we have to 
compute the square root in matrix sense of the diffusion 
matrix D in (^2|). This is not trivial because this is a 
5x5 matrix (rather, a 4 x 4 because it has the first 
column and first row equal to zero) that as elements has 
continuous matrices with indices rr'. However, one can 
resort to the physics of the problem to get a hint on the 
solution. The hint comes from Eqn. (10) that states 
that the time derivative of conserved variables is in the 
form of the divergence of a flux. This suggests that if we 
look at the noise term ( ^7|) it should have the form of a 
divergence. A little of reflection shows that wc can look 
for a noise term of the form 



(48) 

In these expressions, the white noise terms are given in 
terms of independent increments of the Wienner process 
according to (p9|), this is 



'v^dt 



<"(<) 



Vydt 



(49) 



The independent increments of the Wienner process sat- 
isfy 

dWr{t)=dW^^'{t) 
dWr{t)dWy' (f) = [S''^'' S''"'' + S''^''6''''']Srr'dt ift = t' 
dWnt)dW^,' it') = 

dW^^ {t)dW^' (f) = df,f,> 5rr' dt iit = t' (50) 

Whenever correlations of the white noise terms with 
other variables have to be computed, one can use the fact 
that in Ito sense the noise and the variables are uncorre- 
lated at the same time allowing for a decoupling of the 
averages. Then one can use of the following correlations 
for the white noise 

{vntK'''''it')) = [S'"''6'"'' + 6''^''S^"'']6{r - r')6{t - t') 
{€''{tK'.{t'))^Q 

{€{t)4(t'))=5,,,5{v-v')5{t-t') (51) 

Now, we make the ansatz that the matrices B in ( |4^ ) 
are local in space, this is 



B'^^r = B^^^^'Sir - r') 
Kr, - B^^Sir - r') 



In addition, the form of the matrix B must be simple 
because, roughly speaking, its square must be like the 
matrix G in (|4^) which contains only products of Kro- 
necker deltas. Therefore we make the further assumption 
that the matrices B"^'^'^ , B^^ have the form 

B^^''" = a{r)d^pd^,, + b{r) [6a,^6p, + Sa^Jf^p] 
Bf = c{v)5^p (53) 

where a(r), 6(r), c(r) are functions to be specified later. 
These expressions have the correct symmetry properties 
for the matrices B's. In the following, wc will prove that 
this ansatz is correct. 

If we introduce the random stress tensor a'^^{t) and 
random energy flux T^{t) in such a way that 



frit) = 9r„r,"(t) 



(54) 



then we see from (| 
the form 



), (p2), and (pSh that they must have 



afit) = a(r)5„/3tr[77r] + 2b{r)i^f 



(55) 



In order to determine the functions a(r), 6(r), c(r) (and 
check that our ansatz is correct) we compute the corre- 
lation of the noises using the expressions ( ^ ) with the 
definitions ( |55| ) and the properties (pT]). We obtain 

X [5c.p^l5{v-v'){Ab\v)) 
+ dr^dr,S{r - r')(462(r) + 8a(r)&(r) 
+ 6a'{r))] 
{f^{t)fAt'))=-5{t-t') 

X [V2j(r-r')(vj?4fe2(r)) 

+ dr^dr,5{-r - r')(v^'(4fo2(r) + 8a(r)6(r) 
+ 6a2(r)))] 
(/r"(i)/r'(^')> 

-6{t-t') [V25(r-r')(c2(r)) 
V2,5(r-r')(v2462(r)) 
5r„9r,<5(r-r')(v>^(462(r) 
8a(r)6(r)+6a2(r)))l 



{mm')) 

{fHt)fUt')) 



(52) 



These expressions should coincide with 
be rewritten as 

(/r"W/r^(0>='5(t-t')(^rr'> 

{fHt)fUt'))^Sit-t'){Dtt,) 

where the diffusion matrix is given in 
tution of (E^) into (|5^) one obtains 



(56) 
which can 



(57) 

1). After substi- 



8 



{frit)fr^'it')) = ~Ht - t') [5^pVl5{v - r')(2T(r)r;(r)) 
+ ar„9r,'5(r-r')(2T(r) (ry^/S + C(r)))] 

{f^{t)ft,{t'))^~5{t-t') [V?<5(r-r')(v?2r(r),7(r)> 
+ 5r„5r/(r-r') 
X (v^(2T(r)(^(r)/3 + C(r))))] 

(/rW/r'(i')>-(/"W/r'(i')> 

(/rW/r'(i')> --'^(i-i') [V?<5(r-r')(2r(r)2..(r)) 

+ V?5(r-r')(v22r(r),7(r)) 
+ ar„ar,<5(r-r') 

^(2r(r)(rKr)/3 + C(r)))>] (58) 



Compatibility between ( pq ) and ([Sg) then requires that 

Ah^{v) = 2T(r)ry(r) 

3a2(r) + 4a(r)6(r) + 2fe2(r) = T(r) i?7(r) + C(r) 

o 

c2(r) = 2T2(r)K(r) (59) 



We observe then that our ansatz for the matrices B is 
correct because this compatibihty condition can be sat- 
isfied. The solution of this system of equations is 



a(r) + ^6(r) = ±^ir(r)C(r) 



5(r) = ±\/ir(r)r?(r) 



c(r) =T(r)V2K(r) 



(60) 



Therefore, the random stress tensor in ( pSj ) will have the 
following form 



df{t)^^2T{v)ri{v) 



± v/3T(r)C(r)r^'iE^r'W 



(61) 



We have decomposed this tensor as a traceless part plus 
a diagonal part. Both components are statistically inde- 
pendent and therefore the sign ± is indifferent. Actually, 
we can compute the correlation of the random stress ten- 
sor and the result is 



X [(2T(r)r/(r)) [S'^'S^'' + S^^'S^"] 
+ (2r(r)(C(r) - ^,y(r)))<5„^V 



(62) 



In the same fashion, we can define a random heat flux 
according to 



(t) - df (t) {t) = r(r) y^Mv)Vr (t) (63) 



This random heat flux is uncorrelated with the random 
stress tensor and its autocorrelation is given by 

(Jr"(t)jf,(t')> - {2T\r)n{r))6^pS{r - r')S{t - t') (64) 



The expressions (|6lD,(|63|) of the random stress tensor and 
random heat flux in terms of the independent increments 
of the Wienner process are the main results of this paper. 

Now, we have all the elements to write the full 
Langevin equations (|30|). After some rearrangement, the 
Langevin equations can be cast in the form 

atp(r,<) = -V-v(r,i)p(r,t) 
9tg(r, t) = -V- [v(r, i)g(r, t) + a(r, t)] 
dteir, t) = -V- (v(r, i)e(r, t) + v(r, t)-a{r, t) + J(r, t)) 

(65) 

The terms vp, vg, ve represent the purely convective part 
of the fluxes. The total stress tensor a and total heat flux 
J are given by a systematic and a random part 



a{r,t) =a{r,t) +d{r,t) 
J(r,t) = J(r,i) -I- J(r,t) 



(66) 



The systematic parts have the form of Newton's and 
Fourier's constitutive equations 



a = pi 
J = -kVT 



2?7(7-^tr[7]l) -Ctr[7]l 



(67) 



where 7=5 [Vv -I- (Vv)^] is the strain rate tensor. The 
random parts are given by (61),(|63|). 



VII. DISCUSSION 

The Langevin equations for non-linear hydrodynamic 
are the same in structure as the usual macroscopic non- 
linear hydrodynamic equations except that they govern 
the instantaneous fluctuating values of the hydrodynamic 
fields and they have, in addition, a random contribution 
to the stress tensor and heat flux. The form in which 
the thermal noise is described here is similar to the Lan- 
dau and Lifshitz theory |^ . The differences arise on the 
stochastic properties of the random parts in these equa- 
tions. The correlations of the random stress tensor and 
heat flux, Eqn. (|62|),(|6^) are to be compared with those 
given by Landau and Lifshitz p^ : 

2 



+ 2Te,(Ce 



:rieq)Sat3S^,i, 



(J? (i) Jf, (t')> = 2r2^Ke,Ja0,5(r - r')S{t ~ t') (68) 
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Here Teq,r]eq,Ceq-: i^eq BiXB the equilibrium values of the 
temperature and the transport coefficients. It is appar- 
ent that the results ( |62[ ) , (|6j) will coincide with (|8|) only 
if the system is at equilibrium and the fluctuations can 
be considered small, in such a way that the decoupling 
approximation 



where e = In a perturbative expansion in e one 

assumes that the solution of (|7^) can be written as a{t) = 



(T(r)ry(r)) = (r(r))(,?(r)) ^Te^r?,, 



(69) 



can be taken. 

Another point of discrepancy is that the random fluxes 
are not Gaussian. What is Gaussian is the Wienner pro- 
cess, but not the product of the Wienner process with the 
temperature and the transport coefficients. This result 
is consistent with that of Ref. where the assumption 
of Gaussian random fluxes lead to unphysical results. In 
particular, it was found that if the Einstein distribution 
is assumed to be the equilibrium solution, then the vis- 
cosities were inversely proportional to the temperature 
and that the heat conductivity was inversely proportional 
to the temperature squared. This dependences are not 
observed in real systems. Conversely, one could obtain 
physically acceptable transport coefficients by relaxing 
the requirement that the Einstein distribution is the equi- 
librium, which is equally unacceptable. A full discussion 
of this point has been made in Ref. |^ . 

The "size of the fluctuations" is what determines the 



validity of the decoupling approximation (69) and we 
comment on this point now. We note that despite of the 
continuum appearance of the Langevin equations ( |65| ) 
there is an implicit dependence of the random parts on 
the coarse graining, i.e. on the volume of the cells. This 
is apparent in the definition ( p9| ) of the noise appear- 
ing in (|6l|) , (|63|) . In principle, different coarse-graining 
procedures (different discretization with distinct typical 
volumes Vr of the cells) will leave the systematic part 
of the equations invariant, but will change the intensity 
of the noise terms. The 1/y^ dependence is consistent 
with the amplitude of the fluctuations encountered in 
equilibrium ensemble theory. Therefore, the question of 
whether the fluctuations are large or small is a matter 
of the level of description used. Fine resolutions of the 
hydrodynamic fields in which the coarse graining cells 
are small will produce large relative fluctuations in the 
hydrodynamic variables. The existence of characteristic 
length scales imposed on the fluid is what dictates the 
level of coarsening required. In the case of colloidal par- 
ticles, it is observed experimentally that submicronic par- 
ticles diffuse whereas the diffusion for larger particles is 
small and negligible. Therefore, in order to "resolve" sub- 
micronic particles one needs a fine level of coarse graining 
such that the fluctuations are relevant. 

The Landau and Lifshitz fluctuatin g h ydrodynamics 
arises as the small noise expansion of the SDE 
obtained in this paper. This can be seen by writing 
schematically the SDE (|^ with (||) as 



ao{t) + eai{t) + e a2{t) + ■ ■ ■. By substitution in ( |70[ ) and 
cancelling order by order one obtains 

dao{t) = K[&a{t)]dt 

dai{t) ^ K'[aQ{t)]ai{t)dt + B[aQ{t)]dW 
da2{t) = K'[&n{t)]&2{t)dt + ^K"[aoit)]ai{t) 



B'[ao 



.{t)]ai(t)dW{t) 



dt 



(71) 



where the primes denote derivatives. The first equation 
for the dominant contribution ao(t) to the solution of 
the SDE constitutes the deterministic non-linear hydro- 
dynamic equations with fluctuations neglected. The next 
equation for the contribution ai{t) is equivalent to lin- 
earize the equations of hydrodynamics around the de- 
terministic solution and add a Gaussian white noise of 
amplitude B[aQ\. At long times when the system has 
relaxed to equilibrium, these equations are the Landau 
and Lifshitz linear equations governing the equilibrium 
hydrodynamic fluctuations p2[ . Higher order contribu- 
tions (52, . . . are usually not considered, although, in prin- 
ciple, they can be computed systematically p3[. 
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da = K[a\dt + eB[a]dW{t) 



(70) 
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